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1. Introduction 

The Lamb shift in hydrogen was discovered in 1947 |JlJ]. Since then, developments in spec- 
troscopy have led to very precise experimental values for the IS Lamb shift [§, §], making 
it the best test of Quantum Electrodynamics for an atom. On the theory side, much effort has been 
put into the calculation of higher order corrections to match the experimental precision (cf . Ref . [Q] 
for a review of the theory of light hydrogen-like atoms). 

In the perturbative calculation it is important to correctly account for the presence of three 
different small parameters. Za describes the binding effects of an electron to a nucleus with charge 
number Z. Self-interactions of the electron lead to additional powers of a, but not of Z. Therefore, 
it is useful to keep a generic Z to distinguish the two different effects, even though our main 
application is hydrogen where Z is equal to one. The third small parameter is the ratio of electron 
and nucleus masses, m/M. The Lamb shift is of order a{Za) A . Corrections also include logarithms 
of Za and m/M. At present, all second-order corrections in a and Za are known, as well as some 
third order ones [|| g]. 

Another important correction is due to the spatial distribution of the nuclear charge. At the 
moment, the theoretical prediction for hydrogen is limited by the uncertainty in the measurement 



of the proton root mean square charge radius. However, new measurements JlOj , |11[ ] of the Lamb 
shift in muonic hydrogen are expected to improve the knowledge of this parameter, soon. The 
advantage of muonic hydrogen is the greater sensitivity to the proton charge distribution due to the 
larger mass, and thus stronger binding, of the muon. 

Here we present a new calculation of the second-order non-recoil corrections of order a 2 (Za) 5 . 
The full result for these corrections was presented first in Ref. [ |l2| ] and improved in Refs. [jl3|, |l4| ]. 
Our result [|l5|] is compatible with the previous ones and has better precision. 



2. Calculation 

Since we are only interested in non-recoil corrections, we consider the scattering of an on- 
shell electron with an on-shell nucleus, where both external momenta have vanishing space-like 
components. Furthermore, we only have to calculate the leading term in the expansion around 
M — > oo. In this particular case we can simplify the interaction with the nucleus when we consider 
the diagram with direct photon exchange together with the diagram with crossed photon exchange, 
as depicted in Fig. |]. In the infinite mass limit, the nucleus propagators become static propagators. 
However, due to the different momentum flow in the two diagrams, the sum of the two becomes 
the difference of two terms which differ only in the sign of the ie prescription of these propagators. 



X - - 

Figure 1: The effective propagator for the interaction with the nucleus. Wavy and double lines denote 
photons and the nucleus, respectively. 
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(a) 



(6) 



(c) 



Figure 2: Sample diagrams. Solid, wavy, and dotted lines denote electrons, photons, and the effective 
propagator, respectively. 



Neglecting the Dirac structure, we have 
1 1 



M- 



(N + k) 2 -M 2 + ie (N — k) 2 — M 2 + ie 2N-k + ie 2N-k- 



-8(k°) 



(2.1) 



where N = (M,0) and k are the nucleus and photon momentum, respectively. The latter is consid- 
ered to scale like the electron mass, which is much smaller than the nucleus mass. Thus, we have to 
compute three-loop electron self-energy diagrams with two photons and one effective propagator, 
which is defined by Fig |j. Sample diagrams are depicted in Fig. ^[ 

Our calculation proceeds as follows. We use QGRAF flloTj to generate the Feynman diagrams, 
and q2 e and exp []l7|, |l8]] to turn them into FORM [|l|] readable code. Finally, we use MAT AD 3 [|2(j|] 
to do the Dirac algebra and express all diagrams in terms of scalar integrals, using self-made rou- 
tines. The next step is the reduction to so-called master integrals using integration-by-parts identi- 



ties pl| , [22| ]. For this we use the program FIRE []25j], which is a Mathematica implementation 
of the so-called Laporta algorithm [[23|, E4j]. 

Sample master integrals are depicted in Fig. [| We used different methods for their evaluation. 
Simpler integrals like the one in Fig. were computed with the Mellin-B antes method 



27Q , using the Mathematica packages MB [ |28| ] and MBresolve Q29[] - For more complicated 
integrals, we used sector decomposition [[?(], as implemented in the program FIESTA [ |32| , |||]. 
The most complicated master integral is the one in Fig. ||(c). Fortunately, this integral is finite, and 
we were able to derive a Feynman-parameter representation which could be integrated with the 



help of the Cuba library [34]. Results for all master integrals are given in Ref. [15]. 




(a) 



('') 



(c) 



Figure 3: Sample master integrals. Solid and dashed lines denote scalar massive and massless propagators, 
respectively. Dotted double lines denote the delta function. 



3. Results 



In order to present our final result, we split the Feynman diagrams into two subsets: diagrams 
with closed electron loops (cf. Fig. 0(a) and (b)) are denoted by the subscript vp, and diagrams 
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without closed electron loops (cf. Fig. ||(c) and (d)) are denoted by the subscript nvp. Our results 
for the corresponding energy shifts are 

8E vp = a2(Z< * )5 (^) 3 m [0.86281422(3)], (3.1) 



nnP \m) 

8E nvp = t¥*L (t) 3 m [-7.72381 (4)], (3.2) 
where pL = mM/(m+M) denotes the reduced mass of the atom, and n is the principal quantum 



number. Results for individual diagrams can be found in Ref. []15[]. 

The best results for the two subsets prior to our calculation have been published in Ref. [35] 
(cf. Ref. [Q] for references of partial results) and Ref. [|l3|], respectively. Our results agree with the 
previous ones within the error bars. However, we improve the precision by two orders of magnitude 
for SE vp , and a little over one order of magnitude for 8E nvp . 

The total result reads 

8E= a 2 (Za) 5 (tL\ 3 m [_ 6,86100(4)], (3.3) 
and the corresponding energy shifts for the IS and the 2S states in hydrogen are 

8E ls = -296.866(2) kHz, (3.4) 
8E 2S = -37. 1082(3) kHz. (3.5) 
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